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We have found an analytic expression for the multivariate 
generating function governing all n-point statistics of random 
multiplicative cascade processes. The variable appropriate for 
this generating function is the logarithm of the energy density, 
In e, rather than e itself. All cumulant statistics become sums 
over derivatives of "branching generating functions" which are 
Laplace transforms of the splitting functions and completely 
determine the cascade process. We show that the branch- 
ing generating function is a generalization of the multifractal 
mass exponents. Two simple models from fully developed 
turbulence illustrate the new formalism. 
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Multifractals have become a popular tool especially in 
fully developed turbulence to analyze the intermittent 
fluctuations occurring in the energy dissipation field jl],^| . 
They have also been applied to characterize the phase 
space structure of strange attractors in chaotic dynam- 
ical systems ||, diffusion limited aggregation ||, high- 
energetic multiparticle dynamics || , geophysics and self- 
organized criticality ||, to name just a few. 

Among the simplest examples of multifractal processes 
are random multiplicative cascade processes in general 
and the multiplicative binomial process |7j|^] in partic- 
ular. Both the multifractal approach and large devia- 
tion theory analyse these processes in terms of one-point 
statistics. However, because different cascade processes 
like, for example, the p-model ^ and the a-modcl ||] 
are indistinguishable in this framework, such one-point 
statistics appear to be too restrictive. In order to discrim- 
inate among these and other suggested models, a general- 
isation of the multifractal approach to n-point statistics 
is thus called for. The analytic multivariate generating 
function presented in this letter is, we believe, the appro- 
priate generalisation. 

Two-point statistics of random multiplicative cascade 
processes were previously discussed in Refs. |n||llj]. An 
approach to calculate the n-point statistics, or in other 
words the spatial correlations to arbitrary order, within 
these models has been presented in Refs. |L2],|l3): there, 
the multivariate generating function of the spatial cor- 
relations was constructed iteratively from a backward 
evolution equation, leading to a recursive derivation of 
spatial correlations. While in the latter approach, mod- 
els like the p- and a-model become clearly distinguish- 
able, the recursive structure of the multivariate generat- 



ing function appears from a mathematical perspective to 
be less than elegant: an analytic solution is much more 
desireable. The realization of this goal is the subject of 
this letter. 

In order to derive the analytic expression for the multi- 
variate generating function of random multiplicative bi- 
nary cascade processes, we first need to find the appro- 
priate variables and to introduce a convenient labelling. 
The cascade prescription is as follows: A given initial in- 
terval of length 1 and unit energy density successively 
splits into 2, 4, 8, . . . subintervals ("bins"). For the lat- 
ter, we employ a binary labeling: after j cascade steps 
(interval splittings), a specific subinterval is character- 
ized by the sequence kik 2 ■ ■ ■ kj with each ki being ei- 
ther or 1. The bin with the label 00 • • • is the left- 
most one, whereas the label 11 • • • 1 belongs to the bin 
on the far right. The energy density ejp—k- belong- 
ing to a j-th generation bin is redistributed 7 nonuni- 
formly onto the two intervals of the next generation: 



(i+i) Jj) 



fci—fcjO — 1k 1 ---k j € k 1 ---k j i 



and e 



(i+i) 

ki---kjl 



Jj+i) M 



The two multipliers q^-.-kjO and qki- ■■kA are drawn from 
a probability density f>(<3% 1 ...fe 3 0j Qki—kji) which we call 
the splitting function. 

After J cascade steps, the energy density belonging to 
bin (fcife • • • kj) is given by the product of multipliers 
over all previous generations 
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For the logarithm of the energy density this relation be- 
comes of course additive; this observation is crucial for 
the subsequent derivation. 

The structure of a specific random multiplicative cas- 
cade model is completely determined by the splitting 
function p(qo,q\) which determines the distribution of 
energy during each split. Hence, p(qo, qi) also completely 
determines the multivariate statistics of a given cascade. 
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The joint probability density p(lneQ 
find at the same time the logarithm of Eg— oo m the bin 
with label (0 • • • 00), the value lne^ 01 in bin (0 • • • 01), 
and lne^ n in bin (1 • • • 11) can therefore be ex- 
pressed fully in terms of the splitting functions at each 
branching: 

p(lne^ ,...,lne^ 1 ) = 
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notation k, = (fci ■ • • kj) in the subscripts, the multivari- 
ate cumulant of order n is found through 

d n K[X^} 
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This probability density can be converted into a multi- 
variate generating function for cumulants in hie. With 
the definition 

K[4\ , A&] = In /exp( £ Xg.. hj J \ 
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The first three multivariate cumulants are 

C Kl = (lne Kl ) c = (lne Kl ), 
Ck ± ,k 2 = (lne/tj lne« 2 )c 

= (lne«i lne K2 ) - (lne Kl )(lne K2 ), 
= (lneK! lne K2 lne K3 ) c 
= (rne Kl lne K2 lne K3 ) 

- ^ (lne Kl lne K2 )(lneK 3 ) 
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These and higher-order cumulants can be found through 
appropriate derivatives of the b.g.f. using If, for ex- 
ample, the same b.g.f. is used for all branchings, and the 
splitting function is symmetric, p(qo>Qi) = p(<Zij9o); the 
two-point cumulant is 



■kj ' 



(4) 



fcj + i ,. . .,kj—0 

rearranging the terms in the exponent of 
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(F^>r the p-model with a splitting function p(qo,qi) = 
|[£(So — (l+a))+£(?o — (1 — a))]£(<?o+?i— 2), the second- 
order cumulant correlation density is shown in Fig. 1. 

We should point out the scope and limitations of the 
solution (||) . Clearly, the generating function in terms of 
hie is applicable to any functional form of the splitting 



, . ,. a . , « ,i , r,/ c u • , i r function or b.g.f. It does not depend on the number of 

and inserting H2M into (Bl), that the 2-fold integral fac- . , , . 

, branches either: trivariate or even higher-variate spht- 



torizes, so that 
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Q[A ,Ai] 



split- 
ting functions can be implemented. Due to the additive 
nature of the b.g.f. 's, the splitting functions can differ 



, Ai_ x ]= ^ ^ Q[^ki—k -i0' Afci— fc _ii]' (^) fr° m generation to generation and even from branch to 



branch. The only (and important) precondition for the 
applicability of Eq. (||) is that the splitting variables q of 
every branching must be independent of all others^. 



In 
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is the "branching generating function" (b.g.f.) for cumu- 
lants, which governs the behavior of the entire cascade. 
Eqs. (@)-(0) represent the long-sought analytic expres- 
sion for multiplicative cascades: because the multivariate 
cumulant generating function K is the sum of all branch- 
ing generating functions Q, one for every branching, each 
b.g.f. can be solved separately and analytically. 

Multivariate cumulants likewise become sums over cu- 
mulants of the individual branching points. Using the 



1 Nelkin and Stolovitzky |l4) argue that the dependence 
of experimental distributions of splitting variables (multiplier 
distributions) on the position of the subinterval J13] suggests 
that the splitting variables are in fact not statistically in- 
dependent; see also Ref. The effects of this deviation 
from statistical independence on the present formalism, and 
to what extent it matters, remain to be investigated. 
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Figure 1: 

Second-order cumulant C k .k' = C( fcl ... fej ) ! ( fc ;. 
p-model, where k = 1 + Ylj—x kj2 J ^ an d 
k' = l+X)/=i J - Parameters have been set to J = 6 



and 



0.4. 



Moments of the J-generation cascade are found di- 
rectly from the moment generating function Z[A ( -' 7 -'] = 
exp(_ftT[A^' 7 ']). From Eq. (|J), it is obvious that Z factor- 
izes into moment generating functions of the individual 
branchings. 

The key input into the analytical expression (Jq) is the 
branching generating function Q[Ao,Ai] of Eq. (R). Its 
properties uniquely fix the spatial correlations of the bi- 
nary random multiplicative cascade model. For the p- 
model, we get 
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The symmetric a-model is similar to the p-model ex- 
cept that it docs not conserve energy in a cascade 
splitting. Given its splitting function, p(qo,qi) = 

I nLo [%k - (! + «)) + - (1 - «))], its b.g.f. then 
reads 



Q[X , Ai] Q _ modcl = i(A + Ai) ln(l - a 2 ) 

l \^ ( l + a 
2 \T~~a 
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which clearly differs from the b.g.f. (|13|) for the p-model. 
Fig. 2 compares the two branching generating functions. 
Note that for Ao = or Ai = the two expressions ( p"3] ) 
and ( |T4| ) become identical; consequently, the one-point 
cumulants 
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of the p- and a-model are also identical. This is the 
reason why, in a multifractal approach, the two models 
look the same asymptotically. To see differences between 
the two models, one must go to the two-point cumulants 
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Within the p-model, the latter are nonzero for even n 
and zero for odd n > 3. Within the a-model, by con- 
trast, all two-point cumulants vanish since its splitting 
function factorizes: p(qo,qi) = p( f ?o)j'( ( ?i ) ■ We hence see 
that the two-point cumulant moments are sensitive to the 
violation of energy conservation in the splitting function. 

For a binary random multiplicative cascade process to 
qualify as a true multifractal process, the splitting func- 
tion needs to conserve energy, i.e. p(qo, qi) — p(qa)S(qo + 
qi — 2); see Ref. Q for a clarification of this point. If 
energy is indeed conserved, it is possible to link the mul- 
tifractal mass exponents r(q) to the b.g.f. Q[Xo, Ai]. Set- 
ting Xi = and, for simplicity, considering the case of a 
symmetric splitting function p(qo, qy) — p(qi, ?o), Eq. (]?]) 
becomes^] 



Q[X , Ar = 0] = In / dq p{qo)qo = H^°) 
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note that for clarity we have written r(Ao) instead of 
the more familiar notation r(q). Hence, for an energy- 
conserving splitting function, its multifractal mass expo- 
nents follow from Q[Xq, Ai]. The reverse need not be true, 
though. 

Also, for the more general case where energy is not 
conserved in the splitting function, the multifractal mass 
exponents r(q) cannot be deduced in a clean fashion and 
contain less information than the b.g.f. This means that, 
for binary multiplicative cascades, Q[Ao,Ai] can be un- 
derstood as the natural generalization of the multifractal 
mass exponents. 

For a one-dimensional cut through the three- 
dimensional energy dissipation field in fully developed 
turbulence, it is most likely that energy is not conserved 
along the cut. Hence, the splitting function will not 
conserve energy and cannot be identified with the scale- 
invariant multiplier distributions fll6| |. Here it would be 
interesting to find a procedure to infer the proper split- 
ting function from data. Given that the experimentally 



2 The univariate version of relation was previously dis- 
cussed by Novikov |l7]] in connection with the statistics of 
generalised multipliers, the so-called breakdown coefficients. 
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measurable cumulants in In e 



(J) 



are n-fold derivatives 



of Q, the latter can in principle be reconstructed from the 
former. With the help of Eq. (0), the b.g.f. can then be 
inverted into the splitting function via a two-dimensional 
inverse Laplace transformation, 



dx dy p(2e~ x , 2e~ y ) e -(Ao+D*-(A 1+ i)2/ 



?[Ao,Ai]-(A + Ai+2) In 2 



(18) 




p-model 




a-model 

Figure 2: 

Branching generating function Q[Ao, Ai] for the p- and a- 
models for a = 0.4, shown for < A < 10; -10 < Ai < 
10. The intersection of Q with the Ao=0 plane is related 
to the multifractal t(Ao) curve through Eq. jl^). Note 
that these curves are identical for the a— and p-models, 
while otherwise the b.g.f.'s differ significantly. 

Of course, when confronted with real turbulence data, 
the proposed inversion will not be this straightforward as 
the problems of homogeneity Jl8| and statistical depen- 
dence of multipliers ]l^-[T^| have to be taken into account. 

Within the above limitations, we envisage many and 
diverse applications of our analytic solution in many 
branches of physics. Besides fully-developed turbulence, 
the case of high-energetic multiparticle branching pro- 
cesses immediately comes to mind. For the latter, the a- 



and p-models have already been used in this context as 
simulation toy models |p||2(J. Implications in this and, 
for example, random multiplicative process calculations 
in large-scale structure formation in the universe pl| re- 
main to be explored. 

ACKNOWLEDGMENTS 

This work was supported in part by the South African 
Foundation for Research Development. PL acknowledges 
support by APART of the Austrian Academy of Sciences. 



Dedicated to the memory of Peter Carruthers, teacher 
and friend. 

[1] C. Meneveau and K.R. Sreenivasan, J. Fluid Mech. 224, 
429 (1991). 

[2] U. Frisch, Turbulence, (Cambridge University Press, 

Cambridge, 1995). 
[3] T. Halsey, M. Jensen, L. Kadanoff, I. Procaccia and B. 

Shraiman, Phys. Rev. A33, 1141 (1986). 
[4] F. Argoul, A. Arneodo, J. Elezgaray, G. Grasseau and R. 

Murenzi, Phys. Rev. A41, 5537 (1990). 
[5] E.A. De Wolf, I.M. Dremin and W. Kittel, Phys. Rep. 

270, 1 (1996). 

[6] P. Garrido, S. Lovejoy and D. Schertzer, Physica A225, 
294 (1996). 

[7] J. Feder, Fractals, (Plenum Press, New York, 1988). 
[8] C. Meneveau and K.R. Sreenivasan, Phys. Rev. Lett. 59, 
1424 (1987). 

[9] D. Schertzer and S. Lovejoy, in Turbulent Shear Flows 4, 
1985, edited by L.J.S. Bradbury, F. Durst, B. Launder, 
F.W. Schmidt and J.H. Whitelaw, (Springer, 1985), pp. 
7-33. 

[10] M.E. Cates and J.M. Deutsch, Phys. Rev. A35, 4907 
(1987). 

[11] C. Meneveau and A.B. Chhabra, Physica A164, 564 
(1990). 

[12] M. Greiner, P. Lipa and P. Carruthers, Phys. Rev. E51, 
1948 (1995). 

[13] M. Greiner, J. Giesemann, P. Lipa and P. Carruthers, Z. 

Phys. C69, 305 (1996). 
[14] M. Nelkin and G. Stolovitzky, Phys. Rev. E54, 5100 

(1996). 

[15] G. Pedrizzetti, E.A. Novikov and A. A. Praskovsky, Phys. 

Rev. E53, 475 (1996). 
[16] K.R. Sreenivasan and G. Stolovitzky, J. Stat. Phys. 78, 

311 (1995). 

[17] E.A. Novikov, Prikl. Mat. Mekh. 35, 266 (1971); Phys. 

Fluids A2, 814 (1990). 
[18] M. Greiner, J. Giesemann and P. Lipa, Phys. Rev. E56, 

4263 (1997). 

[19] A. Bialas and R. Peschanski, Nucl. Phys. B273, 703 
(1986). 

[20] P. Lipa and B. Buschbeck, Phys. Lett. B223, 465 (1989). 
[21] J. Pando, P. Lipa, M. Greiner and L.-Z. Fang, preprint 
AZPH-TH/97-03, Ap. J. (in press). 



4 



